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ABSTRACT 

We investigate variations of the mass transfer rate in cataclysmic variables (CVs) that are induced by 
nova outbursts. The ejection of nova shells leads to a spread of transfer rates in systems with similar 
orbital period. The effect is maximal if the specific angular momentum in the shell is the same as the 
specific orbital angular momentum of the white dwarf. We show analytically that in this case the nova- 
induced widening of the mass transfer rate distribution can be significant if the system, in the absence 
of nova outbursts, is close to mass transfer instability (i.e., within a factor of ~ 1.5 of the critical mass 
ratio) . Hence the effect is negligible below the period gap and for systems with high-mass white dwarfs. 
At orbital periods between about 3 and 6 hrs the width of the mass transfer rate distribution exceeds 
an order of magnitude if the mass accreted on the white dwarf prior to the runaway is larger than a 
few IQ-^Mq. At a given orbital period in this range, systems with the highest transfer rate should on 
average have the largest ratio of donor to white dwarf mass. We show results of population synthesis 
models which confirm and augment the analytic results. 

Subject headings: novae, cataclysmic variables — stars: binaries: close — stars: evolution — stars: mass 
loss — stars: low mass 



1. INTRODUCTION 

A long-standing problem in the theory of cataclysmic 
variable (CV) evolution is the large scatter of observation- 
ally deduced mass transfer rates for systems with simi- 
lar orbital periods (Patterson 1984, Warner 1987, Warner 
1995). This appears to be in conflict with the well founded 
assumption that mass transfer in CVs is driven by or- 
bital angular momentum losses such as gravitational radi- 
ation and magnetic braking (see, e.g., Rappaport, Joss, & 
Webbink 1982 (RJW), Rappaport, Verbunt, & Joss 1983 
(RVJ), Hameury et al. 1988; for a review see, e.g., King 
1988). 

Binary evolution calculations have shown that CVs 
that form with different initial secondary masses converge 
rather quickly to an essentially unique track in, e.g., the 
orbital period - mass transfer rate (Porb — M) space (e.g., 
RVJ, Kolb & Ritter 1992). This track convergence arises 
as the angular momentum loss rate, and hence the mass 
transfer rate itself, depend on the current system parame- 
ters only. Low-mass main-sequence stars that are subject 
to mass loss rapidly settle at a new equilibrium radius R2 
(Stehle et al. 1996) which in turn is a function of the mass 

3/2 

loss timescale. Therefore the orbital period Porb k R2 
of systems at the equilibrium radius is set by the mass 



transfer rate. 

The track convergence implies that any variation of M 
at a fixed Porb should largely reflect the dependence of 
the angular momentum loss rate J on the white dwarf 
mass Ml. Since —M2/M2 oc J/J, cf. eq (11) below, 
and J has no explicit dependence on Mi for mechanisms 
like magnetic braking which are rooted in the secondary 
star, the expected dependence on Mi at a fixed period 
is oc (Ml -I- M2)^/^/Mi. This is too weak in view of the 
observed scatter of at least an order of magnitude. 

Classical nova outbursts represent an obvious mecha- 
nism that could cause deviations of the transfer rate from 
the secular mean value set by the angular momentum loss 
rate. All CVs are expected to experience some form of 
nuclear burning of the accreted hydrogen-rich material on 
the surface of the white dwarf (WD). In most cases this 
should be non-steady and lead to systemic mass loss, e.g. 
in the form of a classical nova outburst. Specifically, in 
the following we assume that the burning is episodic and 
always accompanied by the ejection of the layers that are 
affected by the burning, i.e., the nova shell. The ejection 
occurs on a very short timescale (< months) and effectively 
leads to an instantaneous change of the orbital separation. 
The magnitude and sign of the change are determined by 
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the amount of angular momentum the shell removes from 
the orbit. In most cases the orbital separation increases, 
with a resulting drop of the mass transfer rate, as the 
Roche lobe is moved suddenly outward by a few atmo- 
spheric pressure scale heights from the photosphere. The 
effect is largest if the angular momentum per unit mass in 
the nova shell is minimal. For an axi-symmetric ejection 
this minimum is just the specific orbital angular momen- 
tum of the white dwarf. We refer to such hydrodynamic 
ejections as ballistic, irrespective of the physical nature 
of the ejection mechanism (e.g., nova explosions, optically 
thick winds from the white dwarf, and so forth). 

If there is additional angular momentum transfer from 
the binary orbit to the nova envelope due to dynamical 
friction of the secondary orbiting inside this envelope the 
increase of the orbital separation upon envelope ejection 
is smaller than in the ballistic case. Very strong frictional 
transfer would even lead to a decrease of the separation. 

The approach taken by MacDonald (1986) to model dy- 
namical friction in the nova envelope gave an extremely 
efficient extraction of angular momentum from the orbit. 
The concomitant mean orbital evolution was greatly accel- 
erated and resulted in both very large mean mass transfer 
rates (> lO~^M0yr~^), and a greater spread of the mean 
mass transfer rate at a given orbital period than in the 
absence of nova outbursts. Such a variation of the mean 
transfer destroys the coherence of evolutionary tracks we 
noted above, and is therefore in conflict with the standard 
model for the CV orbital period gap (see e.g. Kolb 1996). 
MacDonald did not consider the effect of individual out- 
bursts on the mass transfer rate distribution. 

In contrast, Schenker, Kolb, & Ritter (1998), following 
Livio, Govarie, & Ritter (1991) to model dynamical fric- 
tion, found a much weaker angular momentum transfer by 
dynamical friction. Schenker et al. largely dismissed the 
effect of nova outbursts on the mass transfer rate distri- 
bution as unimportant, for ballistic ejections as well as for 
cases with plausible friction. Note that for the remainder 
of this paper we use the terms "mass transfer rate distribu- 
tion" and "mass transfer rate spectrum" interchangeably. 
More precisely, they both describe the distribution in val- 
ues of M at a given orbital period, for either the evolution 
of a single CV or an ensemble of CVs, as a function of 
orbital period. 

A closely related issue is the notion that CVs "hiber- 
nate" after nova outbursts, i.e., there has been the sugges- 
tion that post-novae are very faint, for a time comparable 
to the nova recurrence time (Shara et al. 1986; Prialnik 
& Shara 1987, Livio & Shara 1987). The hibernation idea 
was motivated by space density issues, the need for smaller 
mass transfer rates to obtain violent nuclear runaways in 
theoretical nucleosynthesis models, and by observations 
suggesting that the remnant systems of historical novae 
are faint. 

In this work we reconsider the maximum possible vari- 
ation of the mass transfer rate about the secular mean 
rate due to the ballistic ejection of nova shells. In Sec- 
tion 2 we review the methods by which the mass transfer 
rates in CVs arc deduced. Wc conclude that while some 
of the spread in the M's is due to the uncertainties in 
determining the transfer rates, much of the spread is due 
to a real variation among different systems. In Section 3 



we develop a simple analytical description for the time- 
dependent behavior of M between nova outbursts, as well 
for the probability densities for M at a given orbital pe- 
riod (with technical details deferred to the Appendix). In 
this same section we apply the analytic model to find the 
previously ignored optimum system parameter combina- 
tion that maximizes the spread of mass transfer rates at 
a given orbital period in an observed sample of CVs. In 
Section 4 wc use this prescription in population synthesis 
models to show the overall effect on a theoretical CV pop- 
ulation with standard assumptions about their formation 
and evolution. A significant fraction of the CV population 
lies in the optimum parameter range found in Section 3. 
Finally in Section 5 we draw conclusions concerning the 
likely effects of nova explosions on the spread in values of 
M observed in CVs. 

2. DETERMINATION OF m IN CVS 

Observational determinations of the mass accretion rate 

in cataclysmic variables show a large spread of order 10 
up to 100 in M at a given orbital period, being most pro- 
nounced in the period range of 3 to 6 hr (see Fig. 9.8 in 
Warner 1995). Before presenting our model which may 
partially account for this substantial spread we briefly re- 
view the observational methods used to determine M and 
examine what "spread" may be introduced simply by un- 
certainties in the measurements or the assumptions that 
go into the analysis and interpretation of the data. 

Smak (1989) computed M — My relationships for 
standard thin accretion disks using the usual multi- 
temperature black-body approximations. These relations 
are still often used for determinations of M (see, e.g., 
Robinson et al. 1999). Smak (1994) later added a finite 
disk thickness to account for a bright disk rim and a less 
oblique view of the far side of the disk. He found that 
My increased by only 0.5 magnitudes in most cases for 
moderate to low binary inclinations i, while at large val- 
ues of i and/or high transfer rates, changes in My of up to 
2.5 magnitudes are possible. Smak's results show that for 
high mass accretion rates (10~^ — 10~^ Mq yr~^) a per- 
fectly known value of My corresponds to an uncertainty 
in accretion rate of a factor of ~ 2 — 10 at a given orbital 
period. 

To make use of the M — My relations, the value of My 
for the disk is needed. Observational determinations of 
My are usually based upon the belief that in the optical 
V band, the fiux received is almost entirely from the ac- 
cretion disk. For CVs with orbital periods in the range 
of 3 — 6 hr, this is a fairly safe assumption. Thus, using 
some method of distance determination, and accounting 
for the binary inclination (see Warner 1995 for the con- 
vention used to convert the observed value into a common 
definition of My), one can then estimate M. To the uncer- 
tainties inherent in the M — My relationships one should 
also add the rather substantial uncertainties in the dis- 
tance determinations which are used to find My. 

An independent method of checking whether these M 
estimates are reasonable, which is not based on the deter- 
mination of My , is to utilize X-ray observations and simple 
accretion disk/boundary layer theory. Here the basic as- 
sumption is that some fraction of the accretion luminosity, 
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GMiM / Ri, will emerge in X-radiation. Thus, if the mass 
Ml and radius i?i of the white dwarf are known to some 
degree from observation or by assumption, measurements 
of Lx can be used to estimate M (see Frank et al. 1992). 
Patterson and Raymond (1985a) provided a comparison 
of observed Lx with a suitable observed rnean of M as a 
function of orbital period, derived from M estimates by 
other methods (e.g. Patterson 1984), and found the fol- 
lowing: Using an optically thick boundary layer model the 
two estimates of M agreed to within a factor of 5 at low M 
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Another powerful method to estimate M is photomet- 
ric eclipse mapping. This uses the standard relation 
r^(i?) ~ SGMiM/Snr^ between disk temperature T and 
disk radius r to allow an estimate of M (Horne 1985). If 
multicolor light curves are available, and the effect of red- 
dening can be estimated, this method is independent of 
distance determinations. Smak (1994) showed that for M 
> W^^Mq yT~^ (orbital periods of greater than ~ 5 hr), 
this technique can yield spurious results. However, for typ- 
ical systems in the 3 — 6 hr range, the agreement between 
the observed and theoretically expected T^{R) profiles is 
generally good, with deviations being a factor of 10 or less. 
Baptista (2000) reviews the current level of knowledge in 
this area. 

Yet another method, using equivalent widths of disk 
emission lines, is based on apparent relationships between 
the strength of the disk emission lines and My from the 
accretion disk. For example, Patterson (1984) used the 
equivalent width of H/3 to present a set of accretion disk 
M values. Smak (1989), however, noted that Patterson 
used a mixture of thin and thick accretion disks at var- 
ious binary inclinations and his assumed M values were 
based on an M — My relation by Tylenda (1981) which, 
in turn, used an unrealistically large accretion disk. Thus, 
the determined M values are likely to be underestimated 
in general by factors of 2 — 6. Using the strength of He 
II emission, an assumed white dwarf mass of 0.7 Mq, and 
a 3 A equivalent width for the line, Patterson and Ray- 
mond (1985b) found agreement between the values of M 
predicted by their method, and the values of Af estimated 
by other methods to within factors of 5 — 10 over a broad 
range of M for a number of CVs. 

We find that, in general, the observational determina- 
tions of M have typical associated uncertainties of factors 
of ^ 3 — 5 and in some cases up to a factor of ~ 10. 
From this information we would conclude that the ob- 
served spread in the M distribution, which is found to 
approach factors of 100, especially for the 3 — 6 hr orbital 
period range, is in large part a real effect. 

In the final analysis, however, perhaps the best evidence 
in favor of a wide spread in the values of M at a given 
Porb , is the coexistence of dwarf novae and nova- like vari- 
ables at the same orbital period. The widely accepted 
explanation for dwarf nova eruptions is provided by ther- 
mal ionization instabilities in their accretion disks (Meyer 
& Meyer-Hofmeister 1981; Cannizzo, Shaffer. & Wheeler 
1988; for a review see e.g. Frank et al. 1992) which, in turn, 
require smaller values of M than are inferred for nova-like 
variables which exhibit no such disk instabilities. 



3. ANALYTIC MODEL 

In this section we develop an analytic model for the mass 
transfer rate spectrum that is induced by the ballistic ejec- 
tion of nova shells. For convenience we defer a portion of 
the derivation to the Appendix. After introducing stan- 
dard definitions describing mass transfer via Roche-lobe 
overflow we simply present the main resulting expressions 
and focus on their discussion. 

3.1. Stationary and nonstationary mass transfer 

For simplicity of notation, let M be the magnitude of 
the mass transfer rate, i.e., 

M = -M2 > 0, (1) 

where M2 is the mass of the Roche-lobe filling secondary. 

In general, the transfer rate is a sensitive function of the 
difference between the stellar (photospheric) radius R and 
the Roche lobe radius Rl of the secondary. In the case of 
CVs a good approximation is 

M = Mp^exp|:^^^|. (2) 

Here H is the photospheric pressure scale height of the 
secondary star, and Mph is the mass transfer rate when 
the Roche lobe is situated at the photospheric radius (see, 
e.g., Rittcr 1988, Kolb & Ritter 1990). Typical values of 
Mph for low-mass main-sequence stars are expected to be 
~ lO-^Moyr-i. 

Over intervals that are short compared to the mass loss 
time scale Im = M2/M, the time derivative of the transfer 
rate (2) is 



dt H \ R Rl 



(3) 



if Mph and H are constant and R — Rl <C i?. 

To make use of (3) we introduce the standard decompo- 
sition of radius derivatives into mass transfer- dependent 
and independent terms (e.g. RVJ, Wcbbink 1985, Ritter 
1996), 

'd\nR\ 

(4) 

th 



R ^ M2 



Rl M2 




(5) 



sys 



Here the stellar mass-radius index C describes the adi- 
abatic response of the donor star to mass loss, while 
{dlnR/dt)th denotes the relative radius change due to 
thermal relaxation. (We neglect any changes of the ra- 
dius due to nuclear evolution), (^l is the Roche- lobe index 
and encompasses changes of the Roche lobe due to mass 
transfer and mass loss, and Jgys are "systemic" orbital an- 
gular momentum losses unrelated to mass loss (e.g., due 
to gravitational radiation and magnetic braking) . J is the 
total orbital angular momentum. 

With (4) and (5) equation (3) can be written as 

dM 
dt 

where we have introduced the relative thickness of the pho- 
tosphere of the donor star, 

e = H/R, (7) 



e M2 t *ev 



(6) 
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(for CVs, e ~ 10 the stability denominator 

^ = c-a, (8) 

usually of order unity, and the evolutionary timescale tev, 



te 



d\nR 

dt 



th 




(9) 



sys 



From (3) we see that mass transfer is stationary (dM/dt 
=0) if 

1 = 1^' (10) 
-ft rti, 

and (6) shows that in stationarity the mass transfer rate 
is 

Ms = (11) 

Hence the transfer rate is of order M2/iev, unless the sys- 
tem is close to instability where D is small (see e.g. RJW, 
RVJ, King & Kolb 1995 for a discussion of mass transfer 
instability). 

3.2. Nova outbursts 

To apply this formalism to CVs that experience repeated 
nova outbursts we consider both the actual evolution be- 
tween outbursts, and a mean evolution. 

We assume that the evolution between outbursts is con- 
servative, i.e., none of the transferred mass is lost from the 
system. 

The outburst itself is treated as an instantaneous ejec- 
tion of the nova shell with mass AM^j . If this is related to 
the accreted mass AMig„ needed to trigger the thermonu- 
clear runaway as 

AMej = 7AMig„ = -7AM2 > 0, (12) 
(7 > 0), then the ratio of maximum to minimum values of 
M during an outburst cycle is 



Mf / 1 4 AMign 

= exp ( - 7 



e ' 3 M,, 



exp 



1 4 AMp: 



e 3 M, 



b 



(13) 



where Mf and Mq are the values of M just before and 
just after a nova explosion, respectively, and Mb is the 
total mass of the binary. Equation (13) assumes that the 
ejected mass carries away the mean specific orbital angular 
momentum of the white dwarf. In this case the Roche lobe 
of the secondary always increases as a result of a nova out- 
burst, leading to a sudden drop of the mass transfer rate. 

We assume that the nova outburst recurrence time ij-ec 
is sufficiently short so that the donor mass M2, the relative 
thickness e = H/R of the photosphere and the evolution- 
ary timescale t^v can be considered as constant between 
outbursts. Then, by direct integration of eq. (6), the 
transfer rate between t = 0, the time immediately after 
the outburst, and t = treo the time immediately before 
the next outburst, is 



M{t) 



Mr. 



(14) 



1 + (M,/Mo - 1) exp(-t/eiev) 

(see, e.g., Schenker et al. 1998). Here Mc is the stationary 
mass transfer rate the system would adopt (for conser- 
vative transfer) in the absence of nova outbursts, except 
when it is dynamically unstable. A strict definition of Mc 
is given in the Appendix, see equation (A12). 



3.3. Mass transfer rate speetrum 

Consider now an ensemble of CVs with similar compo- 
nent masses (Mi, M2) and mean mass transfer rate Mm 
at orbital period P. The mean transfer rate is defined by 

r-trec 

Mrr,.= — I M{t)dt (15) 



Mm=T^ [ 
^rec Jo 



and corresponds to the "fiducial" mean evolution where 

the sequence of discontinuous mass ejections is replaced 
by a continuous, isotropic wind, such that Mwind = 

AMej /tree- 

Then the actual (instantaneous) mass transfer rates in 
the ensemble are between Mq and Mf = M(frec)- The 
probability density distribution n(M) in M space is pro- 
portional to 1/M, so that the normalized mass transfer 
rate spectrum reads 



n(M) = -^-^(M) 

Irec dM 



(/n(M)dM = 1). With d</dM 
by manipulating equation (6): 



dt 
dM 



M, 



(^Mc -M^ M 



(16) 

1/ (dM/dt) we obtain 
(17) 



As can be seen from (13) the full ratio of mass transfer 
rates exceeds an order of magnitude if the ejected mass 
is larger than about 1O~^M0. Yet, in many cases, the 
systems still spend most of their time close to the secular 
mean M^, even if the amplitude is large. The correspond- 
ing mass transfer rate spectrum n(logM) oc Mn(M) has 
a narrow spike at large transfer rates, and a long tail with 
low probability density toward small M. An illustrative 
example of M vs. time during the interval between nova 
events is shown in the upper panel of Fig. 1. The lower 
panel in the same figure shows the corresponding proba- 
bility density in M . The cases studied by Schenker et al. 
(1998) had an even narrower distribution than the example 
shown in Fig. 1. 

To judge the effective "width" of the distribution in 
M— space we consider the half width w to be the following 
ratio: 

Mf 



w : 



M1/2' 



(18) 



where M1/2 = M(trec/2) is the median of the mass transfer 
rate distribution (16), i.e. 

"1/2 fMf 

n{M)dM= / n(M)dM. (19) 

Mo JMi/2 

As a working criterion, we adopt the condition «; > 10 
to signal a "wide distribution" . If we define a "stability" 
ratio 

«=|^ (20) 
and a mass accumulation factor 



X = exp 



2e M2 



(21) 
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(a; > 1), then equations (14) and (18) and the above defi- 
nitions can be combined to yield the following expression 
for w: 

= ^il±^. (22) 
^ ^ 1 + x ^ ' 

The stability denominators Dc and Dm are given by (8) 
and refer to the conservative evolution between outbursts 
and the time-averaged evolution, respectively. The stellar 
index is the same for both and D^j while the appro- 
priate Roche lobe indices are given in the Appendix. 
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Fig. 2. — Half width w(x) as a function of mass accumulation fac- 
tor for different stability ratios a (as labelled). The dashed curve 
is the limiting case a = 1. The a increment between the curves is 
0.2 if a > 0. 

The shape of w(x) depends on the value of the parame- 
ter a which describes the stability of mass transfer in the 
system. We require that CVs are dynamically stable with 
respect to the mean evolution, i.e., _D„i > 0. As the nova 
shell carries the mean specific orbital angular momentum 
of the WD we also have < Dm, hence a < 2 (see the 
Appendix). The classical stability limit for conservative 
mass transfer {Dc = 0) is at a = 0, while a <0 describes 
systems where mass transfer would be dynamically unsta- 
ble in the absence of nova outbursts. 

As can be seen from Fig. 2 the character of w{x) changes 
at a = 1 (note that w{l) = 1 in all cases): 



Fig. 1. — Upper panel: Typical evolution of the mass transfer rate 
M between nova outbursts (here for M2 = 0.3Mq, Mi = 0.6Mq, 
AAfign = 3 X 10~^Mq). The time is in units of the recurrence time. 
The upper solid line is the asymptotic value Mc, the lower solid 
line the timc-avcrage Mm = 2 X lO~®M0yr~^. Lower panel: The 
corresponding probability density distribution n(logM). 

As C is essentially a function of stellar mass, and (l is 
a function of the mass ratio q, 

M2 



Q = 



Ml' 



(23) 



only, the stability factor a is fixed for a given system. 
In contrast, the mass accumulation factor x depends on 
the ignition mass. Observational estimates of nova ejecta 

masses, though subject to large uncertainties, tend to give 
much larger values than results of detailed model calcu- 
lations of the thermonuclear runaway on the surface of 
hydrogen-accreting white dwarfs (see, e.g., Starrfield 1999 
and references therein). Hence we examine ui as a func- 
tion of X in order to provide possibly new insight into nova 
ejecta masses based on the spread in CV mass transfer 
rates. 



If 2 > a > 1 the function w(x) reaches a maximum 
at some x — xq, and w{x) 1 for x 00. Wc 
have w < w{xq) < 2 in all cases. For the limiting 
case a = 1, w{x) is monotonically increasing and 
approaches the value w = 2 asymptotically. We 
conclude that for a > 1 the mass transfer rate 
spectrum is always narrow. 



• If a < 1, w{x) is monotonically increasing without 
bound. The slope of w{x) increases with decreasing 
a. We have w{x) = a; for a = 0. The condition 
w > 10 can be met in all cases, for any x > .Xcrit, 
but Xcrit is a function of a (shown in Fig. 3). As 
a is in turn a function of mass ratio and adiabatic 
stellar index, the requirement a < 1 is equivalent 
to a condition of the form q > gcrit- A good 
approximation is 



gcrit - 0.45 + 0.4 C + 



1 



(24) 



6 



Nova-induced mass transfer variations 



10 



4 - 



2 - 




-4 



-3 



-2 



-1 



Fig. 3. — The minimal mass accuinulatioii factor Xcrit 

(a) to obtain 

w{x) > 10. 
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Fig. 4. — Lower limit on the nova trigger mass AMign for nova- 
induced widths of the mass transfer rate spectrum in excess of a 
factor 10. AMign is in units of e ~ 1O~*M0 and shown as a func- 
tion of donor mass or orbital period P/hr = 9{M2/Mq), for 
different white dwarf masses Mi, as labelled. 

The critical conditions q > Qcrit, x > Xait signal- 
ing "wide" mass transfer rate distributions becomes more 
transparent when expressed as constraints on physical sys- 
tem parameters. For a given WD mass the secondary mass 
M2 has to be large enough so that the system is not too 
far from mass transfer instability, and the nova outburst 
trigger mass has to be larger than a certain lower limit, 



AMi„ 



> 



Specifically, the former requirement is 



satisfied when q > qait, where ^crit is given by equation 
(24) and is typically within a factor of ~ 1.5 (smaller than) 



the value of q required for mass transfer stability in the 
absence of nova explosions. Both requirements are sum- 
marized in Fig. 4 where we plot AMiim (in units of cMq) 
as a function of M2 for various Mi. In calculating these 
limits we assumed that 7=1, i.e., the mass of the ejected 
nova envelope is equal to the mass accreted since the last 
nova outburst (this affects Dm) and that the secondary's 
adiabatic mass-radius index C is the same as for a ZAMS 
star with mass M2 (shown in Fig. 5). The figure also gives 
a rough estimate of the orbital period the corresponding 
CV would have, P/hr ~ 9{M2/Mq). 



\r 2 




Fig. 5. — The adiabatic stellar mass radius index as a function 
of stellar mass, for Population I ZAMS stars (from Hjellming 1989). 

The main conclusions to be drawn from Fig. 4 are: 

The nova induced width of mass transfer spectra is 
msignificant below the CV period gap (P < 2 hr). These 
short-period systems are not close enough to instabil- 
ity, i.e., they do not conform to the mass ratio crite- 
rion q > 0.5. For the same reason the nova-induced 
width is negligible whenever the white dwarf is massive 
(Ml > IM0), whatever the orbital period. 

Significant nova-induced widths to the mass transfer 
rate spectrum can occur above the period gap, for systems 
with C/0 white dwarfs of mass 0.55 — O.SMq and secon- 
daries with mass close to the stability limit M2 < 0.7Mi. 
A necessary requirement is, however, that nova trigger 
masses are in excess of a few times 1O~^M0. The ac- 
tual value depends on the photospheric scale height H of 
the donor star, and increases with increasing white dwarf 
mass. It is interesting to note that theoretical outburst 
calculations unanimously predict that the trigger mass de- 
creases with increasing white dwarf mass. Hence we might 
expect that in a realistic CV population only those C/0 
white dwarfs with the smallest mass contribute to an ef- 
fective widening. 

Chiefly, significant widths in M can occur above the 
period gap, but should be least pronounced immediately 
above the gap and for long P. 

In the next section we demonstrate the analjd;ic findings 
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in full population synthesis models. 

4. POPULATION SYNTHESIS MODELS 

In the previous sections, we have shown with an ana- 
lytic model how the mass transfer rate in a CV evolves 
with time during the intervals between nova explosions 
(eq. 14). We then further derived analytic expressions for 
the corresponding probability density distribution in M 
(eqs. 16 and 17), and the half width of this distribution 
(eq. 22). These expressions are straightforward to evalu- 
ate for a given binary system in a particular evolutionary 
state. The behavior of the probability density distribution 
is a sensitive function of the binary properties, including 
the mass ratio and the stability ratio (eq. 20) ; this depen- 
dence is illustrated in Figures 2 — 5. In order to compare 
these results with observational data we would like to know 
what the width in M is as a function of orbital period for 
the ensemble of CV systems that populate the Galaxy. 

In order to explore how broad the range in M would be 
for an observed collection of CVs with different evolution- 
ary histories and binary parameters, we have carried out a 
population synthesis study with the effects of nova explo- 
sions incorporated. The basic approach to the population 
synthesis is discussed in detail in Howell, Nelson, & Rappa- 
port (2001, HNR; see also Kolb 1993, and Kolb 1996, 2001 
for reviews). We start with a large sample (~ 3 x 10^) of 
primordial binaries whose properties arc chosen via Monte 
Carlo techniques. The primary mass is chosen from an 
assumed initial mass function (IMF; Miller & Scalo 1979; 
Eggleton 2001). We adopt a distribution of mass ratios 
for primordial binaries, p{q), which is approximately flat 
(see, e.g., Duquennoy & Mayor 1991). The initial orbital 
period, Porb, is chosen from a function which is constant 
in log Porb (see, e.g., Duquennoy & Mayor 1991; Abt & 
Levy 1985). 

The evolution of these systems is followed to see which 
ones undergo a common envelope phase. In such events, 
the envelope of the evolved giant primary engulfs the sec- 
ondary, leading to a spiral-in episode which leaves the sec- 
ondary in a close orbit with a white dwarf. We adopt 
the standard approach of assuming that some fraction, 
OicE ^ 1, of the initial orbital binding energy is deposited 
into the CE as frictional luminosity (see, e.g., Meyer 
& Meyer-Hofmeister 1979; Sandquist, Taam, & Burkert 
2000). Primordial binaries that are too wide will not un- 
dergo any significant mass transfer and will not lead to the 
formation of CV systems - such wide systems are discarded 
in the present study. Successful pre-cataclysmic variable 
systems to emerge from the first part of the population 
synthesis calculations are those which do undergo a com- 
mon envelope phase and yield a close binary consisting of 
a white dwarf and low mass (< 1 Mq) companion. 

The second part of the population synthesis consid- 
ers those white-dwarf main-sequence binaries for which 
systemic angular momentum losses, or a modest amount 
of evolution by the normal companion star, can initiate 
Roche-lobe contact within a Hubble time. At the start 
of mass transfer from the original secondary to the white 
dwarf the mass transfer may be driven by systemic angular 
momentum losses, or if the secondary has a mass which is 
too large, the transfer may proceed on either a thermal or 
dynamical timescale. We do not consider either of these 



latter two cases, although the former most likely leads to 
a "cousin" of CVs, the supersoft X-ray sources (see, e.g., 
Kahabka & van den Heuvel 1997). Tests for stability at the 
onset of mass transfer and all other phases of the CV evo- 
lution are done in the binary evolution code by requiring 
that D as defined in equation (8) is always positive. The 
binary evolution code is the same as that used by HNR 
which, in turn, is a descendant of the code developed by 
RVJ, and is based on a bipolytropic stellar structure. Each 
of the CV systems is evolved in detail through the mass- 
transfer phase (CV phase) until the donor star has been 
reduced to a negligible mass (typically 0.03 Mq). The ex- 
pression for magnetic braking, which is the principal driver 
of mass transfer for systems above the period gap, is taken 
from RVJ and Verbunt & Zwaan (1981) and adopts values 
for RVJ's parameters of 7rvj = 3 and / = 1.0. In a typi- 
cal population synthesis run, some ~ 2 x 10* such binaries 
are evolved through the CV phase. 

We have modified the treatment of mass transfer in the 
bipolytrope code so that we can explicitly follow the move- 
ment of the Roche lobe through the atmosphere of the 
donor star as the system evolves between nova explosions. 
A few illustrative binary evolution runs were made with 
this modified code (these are not used or presented in this 
paper). However, this type of explicit calculation requires 
very short timesteps (sometimes as short as years) , and is 
therefore very computationally intensive for use in a pop- 
ulation synthesis study. Fortunately, the binary system 
and donor star properties do not evolve significantly dur- 
ing a single interval between nova explosions, and therefore 
much of the analytic formalism developed in this paper can 
be simply carried over to the population synthesis calcu- 
lations. In fact Schenker et al. (1998) showed explicitly 
that CV evolution with the appropriate steady mass-loss 
rate approximates, to a high degree, the average evolution 
of a CV with discrete nova mass loss events. Thus, in the 
population synthesis study we utilize a steady mass loss 
rate wherein the ejected matter in nova explosions carries 
away the specific angular momentum of the white dwarf. 

For each time step in the evolution of a particular CV 
we store the amount of time the system spends within a 
given range of M and Porb- However, for each time step, 
in addition to recording the actual computed value of M 
which represents the case of continuous mass loss, we uti- 
lize equations (14) and (17) to tell us how M is distributed 
in the nova-evolution case. We take from the binary evo- 
lution code instantaneous parameters such as the stellar 
masses, atmospheric scale height, average time to accrete 
an ignition mass, adiabatic stellar index, and so forth, i.e., 
all the parameters necessary to evaluate equations (14) 
and (17) . 

The results are shown first in the M — Porb plane, color 
coded according to the logarithm of the probability of find- 
ing a system in a particular part of this parameter space. 
In Figure 6a the results are shown for the case of a steady 
mass loss instead of the more appropriate discrete nova 
mass-loss events (see also HNR). Note that for systems 
above the period gap the half width of the distribution in 
M at any given value of Porb is no greater than a factor of 
^ 2. In Fig. 6b the same results are shown for the case of 
discrete nova events, where the ignition mass is taken to 
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Fig. 6. — Probability density distribution (PDF) of a CV population in the orbital period - mass transfer rate plane. The color-code refers 
to log(PDF). See text for details. Panel A: Mean evolution, no novae. Panel B: Effect of novae taken into account, trigger masses given by 
eq. 25. Panel C: Effect of novae taken into account, trigger mass 1 X 10"'' = const. Panel D: Effect of novae taken into account, trigger 
mass 3 X lO^''Af0 = const. 




Fig. 7. — Slices in orbital period of the PDFs shown in Fig. 6, on a linear scale. Solid curves: trigger mass 3 x 10 *Mq (Fig. 6d). Dashed: 
no nova case (Fig. 6a). Each panel is labelled with the corresponding orbital period in hrs. 
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be 

AMign ^ 4.4 X 10"'' M© X 

V 109cm; \MqJ yio-^ Mqyt-^ j ^ ' 

is the white dwarf radius). This is an approximate 
analytic fit (cf. Kolb 1995) to results from fully time- 
dependent thermonuclear runaway (TNR) calculations for 
cold (old) white dwarfs by Prialnik & Kovetz (1995), which 

wc use for illustrative purposes only. For the WD radius 
we used the approximation by Nauonberg (1972) 

^ = 0.0112I f - f ^ V'T'^- 

(26) 

We also tested other expressions that fit TNR calcula- 
tions reported in Nomoto 1982, and semi-analytic solu- 
tions to the thermal equilibrium behavior of H-rich shells 
given by Fujimoto (1982). The population synthesis re- 
sults are quite similar in all cases. In the results leading to 
Fig. 6b all CVs with He white dwarf accretors have been 
eliminated. This is due to the uncertainties in our un- 
derstanding of the stability and properties of the nuclear 
burning on the surfaces of such low-mass white dwarfs. 
Note that for the remaining systems, most of which would 
have CO white dwarfs, the width of the M distribution 
at a given Porb is somewhat broader than in the case of 
CV evolution without novae (e.g.. Fig. 6a). For systems 
above the period gap the full width at 50 percent proba- 
bility of the M distribution is about a factor of 3, but is 
still not nearly broad enough to explain the observational 
data on CVs. By contrast, for systems below the period 
gap, the width looks quite broad due to the logarithmic 
color coding; however, a more quantitative look (sec also 
Fig. 7) will show that most systems would have values of 
M closely bunched around the values found in the no-nova 
case. The long colored tail which appears toward lower M 
has quite small probability densities (see e.g. Figure 1), 
consistent with the findings of Section 3. 

Finally, in Fig. 6c and Fig. 6d we show the results for 
the somewhat ad hoc assumption that all novae burn and 
eject 1 X 10"'' Mq or 3 x 10""* Mq of matter. Fig. 6c, for 
dMign = 1 X 1O~^M0, shows a rather uniform width in 
M (above the period gap) of about a factor of 3, which 
is comparable with that for the case where liAfjgn is cal- 
culated from equation (25) (see Fig. 6b). In contrast, the 
distributions in Fig. 6d, which were computed for the case 
where dMig^ = 3 x 10"^Mq, show a width in M that is 
nearly a factor of 10, in reasonable agreement with the 
observations. 

The color-coded images provide a qualitative overview 
of the evolutionary properties expected for CVs in the 
M — Porb plane. However, for a more quantitative view 
we take slices in orbital period and plot ordinary linear 
histograms of the distribution in M for each of 16 inter- 
vals of Porb- The results corresponding to the color image 
of Fig. 6d are shown in Fig. 7 as solid curves. The dashed 
curves (narrower distributions) are for the no-nova cases 
and are shown for purpose of comparison. In each panel, 
the interval in Porb is typically 1/4 hr wide, and centered 



on the value indicated in that panel. A perusal of this fig- 
ure confirms the description of the width in M for various 

regions of Porb as discussed in connection with the color 
figures in the Porb ~ M plane, and the analytic expressions 
for the mass transfer rate spectrum in section 3.3. 

We note that for donors with mass > 0.6 Mq the bipoly- 
trope models underestimate the adiabatic stellar index. 
This might cause a slight overestimate of the width of the 
mass transfer rate spectrum at periods longer than about 
5 hrs. 

5. DISCUSSION 

We have investigated the width of the distribution of 
mass transfer rates for CVs with similar orbital period, 
induced by nova outbursts, as a function of orbital period. 
We used a semi-analytical model that describes the time 
evolution of the mass transfer rate between outbursts to 
probe the parameter space. A detailed population synthe- 
sis model demonstrates the overall effect on a population 
of CVs with standard assumptions about formation and 
evolution. 

CVs where mass transfer is close to instability or un- 
stable, in the absence of nova outbursts, but stabilized 
through outbursts, contribute most to an effective widen- 
ing of the M distribution. Hence the effect is most pro- 
nomiced where this instability occurs, i.e., in systems with 
large mass ratio M2/M1. For this reason the effect is in- 
significant below the period gap (despite possibly large ra- 
tios of maximum to minimum values of M after and prior 
to nova explosions), and at long orbital periods where all 
CVs are dynamically stable. In the intermediate period 
regime, chiefly between 3 and 6 hrs, a surprisingly large 
fraction of CVs arc close to the optimum parameter range 
that maximizes the impact of novae on the spread of M. 
The main contribution is from systems with white dwarf 
masses less than O.SMq. 

We find that the spread of mass transfer rates is sig- 
nificant, i.e., covering an order of magnitude, only if the 
mass accreted between two outbursts exceeds a few times 
1O~^M0. This critical value increases with increasing 
white dwarf mass, and is inversely proportional to the rel- 
ative scale height H/R of the donor star's photosphere. 

Our goal was to determine the maximum effect of nova 
outbursts on the spread of M values. Therefore we as- 
sumed in our analysis that there are no frictional orbital 
angular momentum losses due to the orbital motion of the 
secondary inside the expanding nova envelope. In that 
case the mass transfer rate always decreases as a result of 
an outburst. If there is dynamical friction, the decrease of 
M becomes smaller with increasing friction, and eventu- 
ally turns into an increase (see, e.g., Schenker et al. 1998). 
A similarly large amplitude change in M as in the absence 
of friction, but in the opposite direction, requires an im- 
plausibly strong dynamical friction. 

We assumed further that the mass of the ejected nova 
envelope is equal to the mass the white dwarf accretes be- 
tween two outbursts (i.e., 7 = 1). The effect of changing 7 
is mainly felt through the stability of the systems. A larger 
7 decreases Cm) see equation (AlO), and hence widens the 
parameter space available to systems that are dynamically 
unstable in the absence of nova outbursts but stable with 
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outbursts. Hence any 7 > 1 increases the effect of novae on 
the spread of M, while any 7 < 1 decreases the effect, com- 
pared to the case we have investigated in detail. If 7 = 2 
(i.e., the white dwarf ejects twice as much mass as it has 
accreted) the limiting accreted mass required to cause a 
half width of an order of magnitude for CVs with a 0.6 Mq 
WD decreases from 3.0 x 10"^ to 1.5 x 10"^ Mq. At 
the same time the curves in Fig. 4 extend to somewhat 
smaller values of M2. On the other hand, if 7 = 0.8 (i.e., 
the white dwarf retains 20% of the mass accreted previ- 
ously) the limiting mass becomes 3.5 x 1O~^M0. The 
curves in Fig. 4 move up, those for higher Mi more so 
than those for lower Mi. Systems with Mi = Q.5Mq and 
M2 ^ O.5M0 are dynamically unstable even in the pres- 
ence of nova outbursts. 

Another effect that could reduce the impact of nova out- 
bursts on the mass transfer rate spectrum is wind losses 
from the disk prior to outburst, so that the white dwarf 
accretes less mass than is actually transferred via Roche 
lobe overflow. If this mass carries away the specific orbital 
angular momentum of the white dwarf the effect is that 
the Roche lobe index describing the evolution between two 
outbursts is no longer the conservative one, but a smaller 
one, given by equation (AlO), with 7 now referring to the 
fraction of transferred mass that is lost in the wind. The 
net effect is to stabilize the systems, and hence to reduce 
the parameter space available to nova-stabilized systems. 

Finally, any weighting in favor of bright systems will 
make the maximum of the mass transfer rate spectra seen 
in Fig. 7 more pronounced, i.e., the distribution of M will 
be more skewed. The distributions will therefore appear 
narrower. 

Our results are largely independent of the assumed func- 
tional form of the magnetic braking law. The differential 
effect of nova outbursts on the width of the mass trans- 
fer rate spectrum is the same for any orbital braking law. 
This is significant for the relative distribution of dwarf no- 
vae and novalikes as a function of orbital period. If the 
critical rate for stable disk accretion separates these two 
classes, then a suitably chosen orbital braking strength, 
combined with the nova-induced M spread, can reproduce 
the observed relative distribution. 

A particularly striking detail of the population syn- 
thesis simulations is the effect on the mass transfer rate 
spectrum at periods immediately above the gap (Fig. 7). 
The unwidened distribution is narrow and centered on 
M ~ 1.5 x 10"^MQyr-^ In the standard model, the 
width of the period gap essentially fixes this value, and 
provides a calibration for the magnetic braking strength 
(see, e.g. Kolb 1996). As noted already by Shaffer (1992), 
this rate is lower than the critical rate for disk stability, 
hence suggesting that just above the gap most CVs should 
be dwarf novae. The exact opposite is observed (e.g. Ritter 
& Kolb 1998). Significantly, the nova-widened M distribu- 
tion has a pronounced maximum above the secular mean, 
at M ~ 2.5 x lO-^Mgyr-i (see panel for P = 3.15 hrs 
in Fig. 7), thus allowing a fraction of systems to appear 
as novalikes. Note, however, that the nova widening does 
not affect the mean value of the M distribution at any 
given orbital period. Hence in the nova widened distribu- 
tions shown in Fig. 7 there would still be more dwarf novae 
than novalikes immediately above the period gap. This is 



different if there is an observational selection effect in fa- 
vor of brighter systems. If the selection is strong enough 
dwarf novae that are present in the intrinsic population 
could be effectively suppressed in the observed sample. 

A comment on "hibernation" is in order. In the bal- 
listic ejection cases considered here the post nova transfer 
rate can indeed be much smaller — by one or two orders of 
magnitude — than the secular mean rate (see, e.g.. Fig. 1). 
Generally, after an outburst the systems evolve back close 
to the secular mean within a time /„trcc which is a frac- 
tion /„ < 1 of the nova recurrence time <rec- At periods 
where the nova-induced spread of transfer rates is negli- 
gible, e.g., below the CV period gap, /„ < 0.1, while for 
systems close to dynamical instability /„ ~ 0.5 (see, e.g., 
Fig 1). In other words, the majority of CVs spend a sig- 
nificant fraction of their time close to the secular mean. 
This is the converse of the original hibernation model hy- 
pothesis /„ ^ 1. As recurrence times are generally very 
long, 10^ yrs or longer, this is consistent with observations 
of apparently very faint post nova systems for the few ob- 
served historical novae which obviously had at most a few 
centuries to recover from the outburst. A similar conclu- 
sion was reached by Iben, Fujimoto & MacDonald (1995; 
their Appendix B). 

Our findings rely on the assumption that the majority 
of CV secondaries are unevolved, i.e., main-sequence stars 
(not necessarily in thermal equilibrium) where core hydro- 
gen depletion is not yet very advanced. This is consistent 
with the standard model of CV formation, as employed in 
our population synthesis calculations. There is evidence, 
however, that a significant fraction of CVs do have a some- 
what evolved donor star (Baraffe & Kolb 2000; see also 
Kolb 2001). This should be followed up by future studies. 

We note that Spruit & Taam (2001) have discussed the 
fluctuation of M in the context of a hypothetical circumbi- 
nary disk, but had difficulties finding a convincing model 
that would at the same time preserve the orbital period 
gap of CVs. 

Irradiation (Ritter, Zhang & Kolb 2000, and references 
therein) has been considered as an effect that superposes 
mass transfer cycles, on the thermal time of the donor 
star, onto the secular mean without changing the mean. 
Without additional driving mechanisms, such as the cou- 
pling to a consequential angular momentum loss (e.g. 
McCormick & Frank 1998) all indications suggest that 
irradiation-driven cycles disappear for systems with donor 
mass < 0.7Mq. 

We conclude that nova-induced variations of the mass 
transfer rate may account for the observed spread of M 
above the period gap, but only if the trigger masses 
on intermediate- mass C/0 white dwarfs are in excess of 
~ 2 X 1O~^M0. CVs with the highest transfer rates should 
be close to mass transfer instability, i.e., should have a 
fairly large mass ratio. This constitutes an observational 
test of the nova-induced spread of CV transfer rates. 
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APPENDIX 



ANALYTICAL MODEL: FURTHER DETAILS 

This appendix provides details and derivations for the model developed in Section 3. 

The Roche-lobe index 

The stationary mass transfer rate (11) depends on the Roche lobe index (l- We obtain by comparing the definition 
(5) with the time derivative of 

J = Ml M2./^ (Al) 
V Mb 

{a is the orbital separation, Mt, = Mi + M2 the total binary mass), and by noting that Rl/ci is a function of mass ratio 
q = M2/M1 only. The time derivative of the left hand side of (Al) is the sum Jgys + Jej of systemic losses and losses due 
to mass that is ejected from the system. We parametrize the latter by 

•^ej=Z^^M6, (A2) 

so that v measures the specific angular momentum of the mass lost from the binary, in units of the mean specific angular 
momentum of the orbit. We further assume that the mass lost from the system is a fraction 7 > of the transferred mass, 

Mj, = 7M2 or Ml = (7 - 1)M2. (A3) 
This is consistent with the definition of 7 in (12) if (A3) is understood as the net effect over one nova cycle. 
Using the notation 

/32 = -V^ 
RL/a dq 

this finally gives 

a=7{(2.+ l)g + (-/3.-2)^}+Cc, (A5) 

with (c, the Roche lobe index for conservative mass transfer (7 = 0), given by 

_ , M2-Mi M. 

- + Ml 

If g < 1 we have Rl/c ~ 0.46(M2/M)V3 (Paczynski 1971), hence 

1 Ml 



so that 



and Cc simplifies to the standard expression 



Mb 3 Mb Ml 
Mean evolution 

If the nova ejecta carry the specific orbital angular momentum of the WD we have p = M2/M1, and (A8) simplifies to 

4M2 

C.--73^+Ce. ^ ^ ^ ^ (AlO) 

The corresponding long-term mean mass transfer rate introduced by (15), the "isotropic wind-average", is 

= ^ (All) 

with Dm = C ~ Cm7 cf. (11). The equivalent expression with Dc = C ~ Cc-, 

Me = (A12) 

represents the conservative stationary rate the system would adopt in the absence of nova outbursts if Dc > 0. If the 
conservative evolution is dynamically unstable, i.e., Dc < 0, the quantity Mc as defined in (A12) is still the relevant 
quantity to use in the expression for the time evolution of M, equation (14). 

In the derivation presented here we assume that the thermal relaxation term {d\nR/dt)th which builds up in the donor 

star experiencing a discontinuous mass loss history duo to nova outbursts is the same as in the time-averaged, continuous 
evolution. In other words, tev in (All) is the same as in (14). This has been confirmed by detailed integrations of the 
secular evolution using full stellar models. 
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Outburst amplitude 

The outburst amplitude ratio Mf/Mo of the mass transfer rate (M/ is M at t = tree) can be calculated from the ejected 
envelope mass. Prom (2) we have 

M;_ M(pre-nova) fl/A^X |^ ^^^3^ 



Mo M( post-nova) [ e \ Rl 

where the index "out" refers to the nova outburst. Using (5) the total change (Ai?i,/i?j:,)totai of the Roche lobe radius 
over one outburst cycle (outburst, followed by conservative evolution until the next outburst) can be written as 

)total ( )out"^( ) i„ter ^"^ ^2 ( J ) 1^^^ ' ^ ^ 

while 



fii J,.,., " M, ■ V ^ 

Here "inter" refers to the evolution between two outbursts, and —AM2 = AMig„ > is the mass transferred in this phase. 
Hence we have 

fARL\ ..AM2 „ .AM2 

Using (A16) in (A13) we find 

With (AlO) for C,m this simplifies to equation (13) given in Section 3.2. 

The mass transfer rate spectrum 
The integrated form of the mass transfer spectrum (17) is 

W=eWln|»^l (A18) 



\Mo{Mc-M) 



cf. Equation (14). 

To calculate either (A18) or (17) for systems with specified masses and secular mean rate we need to find Mq. This can 
be obtained as follows. From (14) at t = tree we have 



Mo ^ Mq/M/ -exp{-W/etev} 
Mc 1 - exp{-ti.ec/eiev} 

Using (A17) to replace Mo/Mf and noting that 

tree ^ AMjgn/M^ ^ Dm AMjgn 

etev etev e M2 

(cf. (15) and (All)) we finally find 

Mq _ exp {Dc AMjgn/e M2} - 1 
Wc ~ exp {Dm AMign/eMs} - 1 ' 



(A19) 

(A20) 
(A21) 
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